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This report shows the latest results on the study of the method to determine the angle 03 of the unitarity triangle 
using Dalitz plot analysis of D° decay from B^ — * DK^ process in a model-independent way. We concentrate 
on the case with a limited charm data sample, which will be available from the CLEO-c collaboration in the 
nearest future, with the main goal to find the optimal strategy for 03 extraction. We find that the analysis using 
decays of Dqp only cannot provide a completely model-independent measurement in the case of limited data 
sample. The procedure involving binned analysis of — > DK^ and -0(3770) — > (/f^7r+7r~ ) n {K^tt+tt~ )r> 
decays is proposed, that allows to obtain the 03 precision comparable to unbinned model-dependent fit. 



1. Introduction 



The measurement of the angle fa (7) of the unitar- 
ity triangle using Dalitz plot analysis of the D° — > 
KgTT + TT~ decay from B ± — ► DK^ 1 process, intro- 
duced by Giri et al. [l| and Belle collaboration Q and 
successfully implemented by BaBar 0] and Belle [U, 
presently offers the best constraints on this quantity. 
However, this technique is sensitive to the choice of the 
model used to describe the three-body D° decay. Cur- 
rently, this uncertainty is estimated to be ~ 10° and 
due to large statistical error does not affect the preci- 
sion of </>3 measurement. As the amount of B factory 
data increases, though, this uncertainty will become 
a major limitation. Fortunately, a model-independent 
approach exists (see p|), which uses the data of the r- 
charm factory to obtain the missing information about 
the D° decay amplitude. 

In our previous study of the model-independent 
Dalitz analysis technique [j| we have implemented a 
procedure proposed by Giri et al. involving the divi- 
sion of the Dalitz plots into bins, and shown that this 
procedure allows to measure the phase fa with the 
statistical precision only 30-40% worse than in the 
unbinned model-dependent case. We did not attempt 
to optimize the binning and mainly considered a high- 
statistics limit with an aim to estimate the sensitivity 
of the future super-B factory. 

The data useful for model-independent measure- 
ment are presently available from the CLEO-c exper- 
iment @. The integrated luminocity at the ^(3770) 
resonance decaying to DD available for the analysis is 
400 pb _1 . By the end of CLEO-c operation this statis- 
tics will grow up to 750 pb" 1 Q- This corresponds 
to ~ 1000 events where D meson in a CP eigen- 
state decays to K^tt +, k~ , and twice as much events 
of 0(3770) — > D°D° with both D mesons decaying 
to KgTT + n~ . Both of these processes include the in- 
formation useful for a model-independent ^3 measure- 
ment. In this paper, we report on studies of the model- 
independent approach with a limited statistics of both 
-0(3770) and B data, using both Dcp — > KgTr + ir~ and 



(KgTr + TT )r)(Kgir + ir ) B final states. 



2. Model-independent approach 

The density of D° — > KgTr + ir~ Dalitz plot is given 
by the absolute value of the amplitude f B squared: 



Pd = PD(m 2 + ,m 2 _) = |/d(™+,to?.)| 5 



(1) 



In the case of no CP-violation in D decay the density 
of the D° decay p B equals to 

Pd = \Id\ 2 = P D(m 2 _,m 2 + ). (2) 

Then the density of the D decay Dalitz plot from 
B^ — > DK^ process is expressed as 



Pb± =\f D +r B e^ B± ^f D \ 2 = 

Pd + r 2 B p D + 2^/p D p D (x ± c + y±s), 



(3) 



where x± , y± include the value of (f>3 and other related 
quantities, the strong phase S B of the B ± — ► DK ± 
decay, and amplitude ratio r B : 

x± = r B cos(S B ± 3 ); y± = r B sin(tf B ± cj> 3 ). (4) 

The functions c and s are the cosine and sine of the 
strong phase difference ASd between the symmetric 
Dalitz plot points: 

c = cos(5£)(rn_i, m_) — #D(m_, ra 2 + )) = cos A5 B ; 



s = sin(<5rj(77i^_, m_) — 8 B {m 2 _, m 2 + )) = sin AS 



d ■ 



(5) 

The phase difference AS B can be obtained from the 
sample of D mesons in a CP-eigenstate, decaying to 
Kgir + Tr~. The Dalitz plot density of such decay is 



Pcp = |/d ± Sd\ =Pd+Pd± 2^/p D p D c (6) 

(the normalization is arbitrary) . Decays of D mesons 
in CP eigenstate to KgTr + ir~ can be obtained in the 
process, e.g. e + e~ — > ^(3770) — > DD, where the 
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other (tag-side) D meson is reconstructed in the CP 
eigenstate, such as K + K~ or KgU>. 

Another possibility is to use a sample, where both 
D mesons (we denote them as D and D') from the 
^(3770) meson decay into the KgTr + n~ state [9j. 
Since -0(3770) is a vector, two D mesons are produced 
in a P-wave, and the wave function of the two mesons 
is antisymmetric. Then the four-dimensional density 
of two correlated Dalitz plots is 

Pcorr(m+, m 2 _,m'l,m'*) = Ifpfp ~ I'dId^ = 

I — : ( ? ) 

Pdp' d +Pdp' d - 2^p D p D p' D p' D (cc' + ss), 

This decay is sensitive to both c and s for the price of 
having to deal with the four-dimensional phase space. 

In a real experiment, one measures scattered data 
rather than a probability density. To deal with real 
data, the Dalitz plot can be divided into bins. In what 
follows, we show that using the appropriate binning, 
it is possible to reach the statistical sensitivity equiv- 
alent to the model-dependent case. 

3. Binned analysis with D C p data 

The binned approach was proposed by Giri et al. 
[l|. Assume that the Dalitz plot is divided into 2jV 
bins symmetrically to the exchange m_ <-> m^_. The 
bins are denoted by the index i ranging from —Af to 
Af (excluding 0); the exchange m 2 + «-> m 2 _ corresponds 
to the exchange i <-> —i. Then the expected number 
of events in the bins of the Dalitz plot of D decay from 
B ± -> M± is 

(Ni) = h B [Ki + r%K_i + 2y/K~K~(xc t + y Sl )}, (8) 

where Ki is the number of events in the bins in the 
Dalitz plot of the D° in a flavor eigenstate, h B is the 
normalization constant. Coefficients a and Sj, which 
include the information about the cosine and sine of 
the phase difference, are given by 

/ y/PDpD cos(/S.5D{m\,m 2 _))dV 

d = , , (9) 

/ / PD dV J p D dV 

Si is defined similarly with cosine substituted by sine. 
Here T>i is the bin region, over which the integration 
is performed. Note that Cj = c_j, s, = — s_i and 
c? + s 2 < 1 (the equality cf + sf = 1 being satisfied if 
the amplitude is constant across the bin). 

The coefficients Ki arc obtained precisely from a 
very large sample of D° decays in the flavor eigenstate, 
which is accessible at B-factories. The expected num- 
ber of events in the Dalitz plot of Dq p decay equals 
to 

(Mi) = h C p[Ki + K^ + 2yfK~K~Ci], (10) 



and thus can be used to obtain the coefficient q. As 
soon as the Ci and Si coefficients are known, one can 
obtain x and y values (hence, 03 and other related 
quantities) by a maximum likelihood fit using equa- 
tion ©. 

Note that now the quantities of interest x and y 
(and consequently 03 ) have two statistical errors: one 
due to a finite sample of B ± — * DK ± data, and due to 
Dcp — + Kgir + ir~ statistics. We will refer to these er- 
rors as B-statistical and Dcp-statistical, respectively. 

Obtaining Sj is a major problem in this analy- 
sis. If the binning is fine enough, so that both the 
phase difference and the amplitude remain constant 
across the area of each bin, expressions (j9|) reduce to 
d — cos(ASd) and Si = shi(A<5d), so s, can be ob- 
tained as Si = ±yT— cf. Using this equality if the 
amplitude varies will lead to the bias in the x, y fit 
result. Since Ci is obtained directly, and Sj is over- 
estimated by the absolute value, the bias will mainly 
affect y determination, resulting in lower absolute val- 
ues of y. 

Our studies [5( show that the use of equality cf + 
sf = 1 is satisfactory for the number of bins around 
200 or more, which cannot be used with presently 
available Dcp data. It is therefore essential to find 
a relatively coarse binning (the number of bins be- 
ing 10-20) which a) allows to extract Sj from Cj with 
low bias, and b) has the sensitivity to the 03 phase 
comparable to the unbinned model-dependent case. 

Fortunately, both the a) and b) requirements ap- 
pear to be equivalent. To determine the B-statistical 
sensitivity of a certain binning, let's define a quantity 
Q — a ratio of a statistical sensitivity to that in the 
unbinned case. Specifically, Q relates the number of 
standard deviations by which the number of events 
in bins is changed by varying parameters x and y, to 
the number of standard deviations if the Dalitz plot 
is divided into infinitely small regions (the unbinned 
case): 




(11) 

where f B = Id + (x + iy)fp, N = J^. \f B \ 2 dV. 

Since the precision of x and y weakly depends on 
the values of x and y Q, we can take for simplicity 
x = y = 0. In this case one can show that 

Q%= y =o = + s 2 i)Ni /J2 N i (12) 

i I i 

Therefore, the binning which satisfies cf + sf = 1 (i.e. 
the absence of bias if Si is calculated as y/l — cf ) also 
has the same sensitivity as the unbinned approach. 
The factor Q defined this way is not necessarily the 
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best measure of the binning quality (the binning with 
higher Q can be insensitive to either x or y, which 
is impractical from the point of measuring 3 ), but 
it allows an easy calculation and correctly reproduces 
the relative quality for a number of binnings we tried 
in our simulation. 

The choice of the optimal binning naturally depends 
on the D° model. In our studies we use the two- 
body amplitude obtained in the latest Belle 03 Dalitz 
analysis [J|. 

From the consideration above it is clear that a 
good approximation to the optimal binning is the 
one obtained from the uniform division of the strong 
phase difference ASp,- In the half of the Dalitz plot 
m 2 ^ < m 2 _ {i.e. the bin index i > 0) the bin T>i is 
defined by condition 

2tt(z- l/2)/Af< AS D (m%,m 2 _) < 2w(i + 1/2) /Af, 

(13) 

in the remaining part (i < 0) the bins are defined 
symmetrically. We will refer to this binning as ASp- 
binning. As an example, such a binning with Af = 8 is 
shown in Fig. Q] (left). Although the phase difference 
variation across the bin is small by definition, the ab- 
solute value of the amplitude can vary significantly, so 
the condition c? + sf = 1 is not satisfied exactly. The 
values of Cj and Sj in this binning are shown in Fig. [3] 
(bottom left) with crosses. 

Figure [T] (right) shows the division with Af = 8 
obtained by continuous variation of the Afe-binning 
to maximize the factor Q. The sensitivity factor Q 
increases to 0.89 compared to 0.79 for Afo-binning. 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

m; (GevV) m; (GevV) 



Figure 1: Divisions of the D a — ► Kgiv + iT~ Dalitz plot. 
Uniform binning of ASd strong phase difference with N = 
8 (left), and the binning obtained by variation of the latter 
to maximize the sensitivity factor Q (right). 

We perform a toy MC simulation to study the sta- 
tistical sensitivity of the different binning options. 
We use the amplitude from the Belle analysis [4] 
to generate decays of flavor D°, Dcp, and D from 
B ± -> DK ± decay to the K° s 7r + 7r final state accord- 
ing to the probability density given by (5), © and 
©, respectively. To obtain the B-statistical error we 
use a large number of D° and Dcp decays, while the 
generated number of D decays from the B ± — > DK ± 
process ranges from 100 to 100000. For each number 
of B decay events, 100 samples are generated, and 



the statistical errors of x and y are obtained from the 
spread of the fitted values. A study of the error due 
to Dcp statistics is performed similarly, with a large 
number of B decays, and the statistics of Dqp decays 
varied. Both errors are checked to satisfy the square 
root scaling. 

The binning options used are AJu-binning with 
Af = 8 and Af = 20, as well as "optimal" binnings with 
maximized Q obtained from these two with a smooth 
variation of the bin shape. Note that the "optimal" 
binning with Af — 20 offers the B-statistical sensitiv- 
ity only 4% worse than an unbinned technique. For 
comparison, we use the binnings with the uniform di- 
vision into rectangular bins (with Af = 8 and Af = 19 
in the allowed phase space, the ones which are denoted 
as 3x3 and 5x5 in Q). 

The B- and Upp-statistical precision of different 
binning options, recalculated to 1000 events of both 
B and Dcp samples, as well as their calculated values 
of the factor Q, are shown in Table fl] In the present 
study we use the errors of parameters x and y rather 
than 03 as a measure of the statistical power since 
they are nearly independent of the actual values of 
03, strong phase 5 and amplitude ratio rp- The er- 
ror of 03 can be obtained from these numbers given 
the value of rp- The factor Q reproduces the ratio of 

the values yjl/cr 2 + 1/cr 2 , for the binned and unbinned 

approaches with the precision of 1-2%. While the bin- 
ning with maximized Q offers better i?-statistical sen- 
sitivity, the best Dcp-statistical precision of the op- 
tions we have studied is reached for the A^-binning. 
However, for the expected amount of experimental 
data of B and Dqp decays the B-statistical error dom- 
inates, therefore, slightly worse precision due to Dqp 
statistics does not affect significantly the total preci- 
sion. 

We have considered the choice of the optimal bin- 
ning only from the point of statistical power. How- 
ever, the conditions to satisfy low model dependence 
are quite different. Since the bins in the binning op- 
tions we have considered are sufficiently large, the re- 
quirement that the phase does not change over the 
bin area is a strong model assumption. We have per- 
formed toy MC simulation to study the model depen- 
dence. While the binning was kept the same as in 
the statistical power study (based on the phase dif- 
ference from the default D° amplitude), the ampli- 
tude used to generate D°, Dqp and B ± — > DK ± 
decays was altered in the same way as in the Belle 
studyof the model-dependence in the unbinned anal- 
ysis [J]. As a result, the same bias of A03 ~ 10° is 
observed as in unbinned analysis. The bias in x and y 
if demonstrated in Fig. [5) We remind that the cause 
of this bias is a fixed relation between the Ci and Si. 
Therefore, proposed binning options, although provid- 
ing good statistical precision, arc not flexible enough 
to provide also a low model dependence. To minimizie 
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Table I Statistical precision of (x,y) determination using different binnings and with an unbinned approach. The errors 
correspond to 1000 events in both the B and Dcp ((Kstt + it~) 2 ) samples. 



Binning 


Q 


B-stat. err. 


Dcp-stat. err. 




7r ) 2 -stat. err. 




a y 


Cx 


a y 


&x 


a y 


Af = 8 (uniform) 


0.57 


0.0331 


0.0600 


0.0053 


0.0097 


0.0145 


0.0322 


A/" = 8 (AS D ) 


0.79 


0.0273 


0.0370 


0.0042 


0.0072 


0.0050 


0.0095 


Af =8 (optimal) 


0.89 


0.0232 


0.0324 


0.0058 


0.0114 


0.0082 


0.0114 


Af = 19 (uniform) 


0.69 


0.0274 


0.0549 


0.0042 


0.0112 






Af = 20 (A5d) 


0.82 


0.0266 


0.0350 


0.0048 


0.0074 






Af = 20 (optimal) 


0.96 


0.0223 


0.0290 


0.0078 


0.0110 






Unbinned 




0.0213 


0.0279 












-0.04 -0.02 0.02 0.04 -0.04 -0.02 0.02 0.04 
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Figure 2: Toy MC study of the analysis using Dcp data. 
Difference between the fitted and generated x (left) and 
y (right) values. Result of the toy MC study with ASd 
binning, 10 5 B decays and 10 4 Dcp decays. Histogram 
shows the fit result with the same D° decay amplitude 
used for event generation and binning, the points with the 
errors bars show the case with different amplitudes. 



the model dependence, the bin size should be kept as 
small as possible, therefore, uniform binning is more 
preferred. 

In a real analysis, one can control the model error 
by testing if the amplitude used to define binning is 
compatible with the observed Dqp data. This can 
be done, e.g., by dividing each bin and comparing 
calculated values of Ci in its parts, or by comparing 
the expected and observed numbers of events in each 
bin. The first results by the CLEO-c collaboration are 
available Q that show the good agreement of experi- 
mental data with Ci calculated from two-body ampli- 
tude for A<5£>-binning. 

We conclude that the method of </>3 determination 
using only Dcp data is only asymptotically model- 
independent, since for any finite bin size the calcu- 
lation of Si is done using model assumptions of the 
ASd variations across the bin. Increasing the Dcp 
data set, however, allows to apply a finer binning and 
therefore reduce the model error due to the variation 
of the phase difference. 



4. Binned analysis with correlated 

D° -> K° s 7i + n- data 

The use of the ^(3770) decays where both neutral 
D mesons decay to the Kgir + Tr~~ state allows to signif- 
icantly increase the amount of data useful to extract 
phase information in D° decay. It is also possible to 
detect events of V>(3770) -> {K s 7r + n-) D (K L 7r+TT-) D , 
where is not reconstructed, and its momentum 
is obtained from kinematic constraints. The num- 
ber of these events is approximately twice that of 
(KgTr + tt~ ) 2 due to combinatorics. However, it is im- 
possible to simply combine these samples since the 
phases of the doubly Cabibbo-suppressed components 
in D° -> if°7T+7r- and TP -> K° L n + n- amplitudes 
are opposite [§]. In the analysis of B data only 
K s n + n~ state can be used, but it is possible to utilize 
K l tt + it~ data to better constrain the D a — > A'g7r + 7r~ 
amplitude using model assumptions based on SU(3) 
symmetry. In what follows, we will consider the use 
of Kgir + iT~ data only. 

In the case of a binned analysis, the number of 
events in the region of the (KgTr + ir~) 2 phase space 
is 

(M)y = h co „[KiK-j + K-iKj- 

(14) 

2^/KiK-iKjK-j(acj + SiSj)]. 

Here two indices correspond to two D mesons from 
^(3770) decay. It is logical to use the same binning as 
in the case of Dqp statistics to improve the precision 
of the determination of Cj coefficients, and to obtain 
Si from data without model assumptions, contrary to 
Dcp case. The obvious advantage of this approach is 
its being unbiased for any finite (Kgir + ir~) 2 statistics 
(not asymptotically as in the case of Dq p data) . 

Note that in contrast to Dcp analysis, where the 
sign of Si in each bin was undetermined and has to be 
fixed using model assumptions, (Kg tt + tt~) 2 analysis 
has only a four-fold ambiguity: change of the sign of 
all Ci or all s,-. In combination with Dqp analysis, 
where the sign of Cj is fixed, this ambiguity reduces 
to only two-fold. One of the two solutions can be 
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chosen based on a weak model assumption (incorrect 
Si sign corresponds to complex-conjugated D decay 
amplitude, which violates causality requirement when 
parameterized with the Breit-Wigner amplitudes) . 

Coefficients Cj, Sj can be obtained by minimizing 
the negative logarithmic likelihood function 

- 2 log C = -2 log P(My, (M)y), (15) 

where P(M, (M)) is the Poisson probability to get M 
events with the expected number of (M) events. 

The number of bins in the 4-dimcnsional phase 
space is 4JV 2 rather than 2Af in the Dqp case. Since 
the expected number of events in correlated iT<j7r + 7r~ 
data is of the same order as for Dqp, the bins will be 
much less populated. This, however, does not affect 
the precision of Cj, Si determination since each of the 
free parameters is constrained by many bins. 

The coefficients Ci, s,; obtained this way can then 
be used to constrain x, y with the maximum likeli- 
hood fit of the B decay data using Eq. [5J To correctly 
account for the errors of Cj, s% determination, this like- 
lihood should include distributions of these quantities, 
in addition to Poisson fluctuations in B data bins. A 
more convenient way is to use the common likelihood 
function, covering both B and (KgTT + ir~ ) 2 data: 

-2 log C = -2 log P(M l3 , (M)y-)- 

lJ (16) 
2^1ogP(iV l ,(iV) J ), 

i 

with x, y, hp, h corr , Ci and Si as free parameters. This 
approach is also more optimal in the case of large B 
data sample, since it imposes additional constraints 
on Ci, St values. 

The toy MC simulation was performed to study 
the procedure described above. Using the amplitude 
from the Belle analysis, we generate a large number of 
D — > Kgir + Tr~ decays and several sets of (.ftr<j7r + 7r~) 2 
decays (according to the probability density given by 
0) and B decays © . We use the same binning op- 
tions as in Dcp study with M = 8. The combined 
negative likelihood (|16[) is minimized in the fit to each 
toy MC sample. We constrain c 2 + sf < 1 in the fit to 
avoid entering unphysical region with negative num- 
ber of events in the bin. The number of (.ftr<j7T + 7r~) 2 
and B decays ranges from 10 3 to 10 5 . The errors of 
x and y parameters are calulated from the spread of 
the fitted values. If the number of (Kgir + ir~) 2 decays 
is comparable or larger than the number of B decays, 
the x and y errors can be represented as quadratic 
sums of two errors, each scaled as a square root of 
(KgTr + ir~) 2 and B statistics, respectively. However if 
the number of B decays is large, the errors of Ci and 
Si depend also on B decay statistics, so separating the 
total error into B- and (i4T^7r + 7r _ ) 2 -statistical errors 
becomes impossible. 



200 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 



Figure 3: Toy MC study of the analysis using (7fg7r + 7r - ) 2 
data. Top plots: difference between the fitted and gen- 
erated x (left) and y (right) values. Result of the toy 
MC study with A8n binning, 10 5 B decays and 10 4 
(Ksiy + it~) 2 decays. Histogram shows the fit result with 
the same D° decay amplitude used for event generation 
and binning, the points with the errors bars show the case 
with different amplitudes. Bottom plots: coefficients d, 
Si obtained in the fit to toy MC sample. Different colors 
correspond to different bins. Cases with the same ampli- 
tude (left) and different amplitudes (right) used for event 
generation and binning. 



The best (iT!y7r + 7r~) 2 -statistical error is obtained 
for Afe-binning and recalculated to 1000 events yields 
o~ x — 0.0050, a y = 0.0095, which is only slightly worse 
than the error obtained with the same amount of Dcp 
data (see Table U for comparison). We also check that 
significant change of the model used to define the bin- 
ning does not lead to the systematic bias (although 
it does decrease the statistical precision). Figure |3] 
demonstrates the precision of the determination of Cj, 
Si coefficients in our toy MC study and the absence of 
the systematic bias for both x and y when the model 
is varied. 

The numbers of (Kgir + ir~) 2 and Dcp decays in 
^(3770) data are comparable, so are the statistical 
errors due to ^(3770) data sample for the two ap- 
proaches. The same binning can be used in both ap- 
proaches, therefore improving the accuracy of Cj deter- 
mination. The approach based on (Kgir + ir~) 2 data 
allows to extract both c, and without additional 
model uncertainties, so it can be used to check the 
validity of the constraint cf + s 2 = 1 and therefore to 
test the sensitivity of the particular binning. 
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5. Conclusion 

We have studied the model-independent approach 
to 3 measurement using B ± — ► DK ± decays with 
neutral D decaying to JsT<j7r + 7r~. The analysis of 
-0(3770) — > DD data allows to extract the informa- 
tion about the strong phase in D° — > ATg7r + 7r~ de- 
cay that is fixed by model assumptions in a model- 
dependent technique. We specially consider the case 
with a limited -0(3770) — > data sample which will 
be available from CLEO-c in the nearest future. 

In the binned analysis, we propose a way to obtain 
the binning that offers an optimal statistical precision 
(close to the precision of an unbinned approach) . Two 
different strategies of the binned analysis are consid- 
ered: using Dqp — > Kgir + ir~ data sample, and using 
decays of 0(3770) to (K%ir + ir-) D (K%n + ir-) D . The 
strategy using Dqp decays alone cannot offer a com- 
pletely model-independent measurement: it provides 
only the information about Ci coefficients, while Si 
for low Dcp statistics has to be fixed using model 
assumptions. However, as the Dqp data sample in- 
creases, model-independence can be reached by reduc- 
ing the bin size. The strategy using the 0(3770) — > 
(K n + 'K~)E>{K®Tt + Tr~)Ei sample, in contrast, allows 
to obtain both Cj and Sj with an accuracy compara- 
ble to Dcp approach. Both strategies can use the 
same binning of the D° — > Kgir + Tr~ Dalitz plot and 
therefore can be used in combination to improve the 
accuracy due to 0(3770) statistics. 

The expected sensitivity is obtained based on the 
D° decay model from Belle analysis. For the CLEO-c 
statistics of 750 pb" 1 (~ 1000 D C p and (K^tt+tt-) 2 
events) the expected errors of parameters x and y due 
to 0(3770) statistics are of the order of 0.01. For rs = 
0.1 it gives the 03 precision — o~ Xt y/{-\/2rB) — 5°, 
which is far below the expected error due to present- 
day B data sample. Further improvement of 03 pre- 
cision will require larger charm dataset, which can be 
provided by BES-III experiment [Io| . 



In our study, we did not consider the experimental 
systematic uncertainties e.g. due to imperfect knowl- 
edge of the detection efficiency or background com- 
position. We believe these issues can be addressed 
in a similar manner as in already completed model- 
dependent analyses. 
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